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Recently, it was found in QED that the generation of a dynamical electron mass in a strong 
magnetic field is significantly enhanced by the perturbative electron mass. In the present paper, the 
related question of a possible enhancement of the dynamical quark mass in an external magnetic 
field and with a bare mass term is investigated in the Nambu-Jona-Lasinio model. 



In the recent paper [l[ some aspects of the well-known magnetic catalysis effect were studied in quantum electro- 
dynamics (QED). In particular, it was shown that an enormously high external magnetic field B me of the order of 
10 82 G would be needed to create dynamically the common experimental value m e w 0.5 MeV of the electron mass 
m e in the massless QED. At the same time, if the bare mass of the theory is nonvanishing and corresponds to the 
experimental value m e , then in the presence of an external magnetic field of the same value B mi! the dynamical mass 
of an electron is enhanced to the value almost ten times larger than m e . (The behaviour of the dynamical electron 
mass in a strong magnetic field in massive QED was also considered previously 0].) As it was claimed in such a 
significant enhancement of the dynamical electron mass in QED is a new effect that can find applications in astro- 
physics and cosmology, and it deserves to be investigated in more detail, and especially also in QCD. In particular, 
as it was predicted in [jj, in the case with much smaller and realistic magnetic field values around 10 15 G, the typical 
magnetic fields of compact stars, a few percent measurable increase in the dynamical electron mass still exists. 

It has been known that the magnetic catalysis effect, i.e. the spontaneous breaking of the chiral symmetry induced 
by an external magnetic field B, is an universal phenomenon, which takes place in different physical models (see, e.g., 
the reviews Q as well as the original papers IS B fl S (IS and references therein). Thus, the following natural 
question arises: Does it mean that the enhancement effect has an universal character as well? In this 
paper, we study this problem in the framework of the Nambu - Jona-Lasinio (NJL) model with two quark flavors. 

In four-dimensional spacetime and at B = the system is described by the following Lagrangian: 



C = q{^d v - m }q + G[{qqf + {qi^rqf 



(1) 



where the quark field q = qi a is a flavor doublet (i = 1, 2 or i = u, d) and a color triplet (a = 1, 2, 3 or a — r, g, b) 
as well as a four-component Dirac spinor; r a stands for the Pauli matrices. It is supposed here that up and down 
quarks have an equal current (bare) mass mo- Clearly, at mo = this Lagrangian is invariant under the continuous 
chiral SU(2)ixSU(2)^ group as well as under the discrete chiral transformation, q — > ij 5 q. At the tree level, the 
Lagrangian (TT]) contains two free model parameters, the coupling constant G and the bare quark mass toq. However, 
when including quantum effects (quark loops), one should regularize the corresponding loop integrals, for example, by 
cutting off the three-dimensional momentum space, i.e. supposing that \p\ < A. Thus an additional free parameter, 
the cutoff A, appeares in the model. In the mean field approximation the effective potential of the model ((T|) looks 
like (see, e.g., Gllll) 
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where m is the dynamical quark mass, which is connected with the bare mass mo and the vacuum expectation value 
of quark fields (qq) through the relation 

m = mo — 2G(qq). 

Note that it depends on the model parameters G, mo, A and is determined by the gap equation 

m — mo 3m 
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Evidently [TTI . [l2j], at m = the dynamical quark mass m is a nonzero quantity only at G > G cr n = 7r 2 /(6A 2 ) 
(in this case the chiral symmetry of the model is spontaneously broken down). However, it follows from © that at 
mo = and G < G cr n we have m = 0, and the chiral symmetry remains intact in this case. If mo 7^ 0, then ra^O 
for arbitrary values of G. Below, one can find some values of m vs mo in the second line of Tables I, II for G < G cr it- 
The influence of an external constant and homogeneous magnetic field B on the properties of the NJL-type models 
was already considered in refs. To obtain the corresponding Lagrangian, it is necessary to perform 



2 



in |T|) the following replacement: d v — > d v + iQA„, where A v is a vector-potential of an external magnetic field B, 
and Q — diag(ei,e 2 ) is the electric charge matrix of quarks. Here e\ — 2|e|/3 and e 2 = — | e| /3 (e is the electric 
charge of electrons) are the electric charges of u- and d- quarks, respectively. At too = 0, the resulting Lagrangian 
is still symmetric with respect to the discrete symmetry q — > ij 5 q, but it is no more invariant under the continuous 
chiral symmetry SU(2)ixSU(2)^j because of the difference in quark electric charges. Clearly, at B ^ the effective 
potential of the model is also modified (see, e.g., [TH) and looks like 

V eS (M; B) = V(M) J2 ^T^'t- 1 ' **) - 5^ - *<] lnx * + f }' (4) 

i=l 

where Xi = M 2 /(2\eiB\) for each i = 1,2, £'(— l,a;)= a;)/dz/| v= _i (C(^, %) is the generalized Riemann zeta- 

function [IB]), and y(M) is the effective potential ([2]) with m replaced by M. The quantity M = M(too, B) in Q is 
the dynamical quark mass which is the solution of the gap equation 

— V eS (M; B) = m V{M) - h{M) - I 2 (M) = 0, (5) 

where 

h{M) = iM } e f\ \nT{ Xl ) - iln(27r) + x t - \{2 Xl - l)hx Xi } (i = 1,2) (6) 

and r(x) is the Euler gamma-function [l(|. In what follows, we suppose that £? is a nonnegative quantity, and take 
into consideration that M(too,0) = to. 

At zero bare mass too = the influence of an external magnetic field on the phase structure of the model |T]) was 
already investigated, e.g., in [l2|, [H, EH EH- In particular, it was shown there that at G < G cr n and B = the 
global minimum of the effective potential © lies at the point to = 0, so that the chiral symmetries, both continuous 
and discrete, are not broken down. However, at arbitrary small values of B and G < G cr n the global minimum of 
the effective potential ((H) of the system is shifted to a nontrivial point. As a result, in this case the spontaneous 
breaking of the discrete chiral symmetry 1 is induced by an external magnetic field B ^ (magnetic catalysis effect). 
Moreover, a dynamical quark mass M = M(mo = 0, B), which is the solution of the equation ([5]) at too = 0, is also 
generated. Note, at G > G cr u and too = the dynamical chiral symmetry breaking in the NJL model takes place 
even at B = due to a rather strong interaction in the quark-antiquark channel. 

Now, we have at our disposal all necessary formulas in order to solve the question raised at the beginning of the 
paper. Clearly, the two different possibilities should be studied, G < G cr i t and G > G cr i t - 

The weak coupling regime (G < G cr u). Since in this case, as in QED, the magnetic catalysis effect takes place in 
the NJL model, i.e. at too = a nonzero dynamical quark mass is induced by an external magnetic field, we are going 
to proceed in the spirit of the paper [lj. For simplicity, let us put A = 1 GeV, i.e. G cr u ~ 1.65 GeV -2 , and consider, 
for illustrations, two values of the coupling constant, G — 0.5 GeV -2 and G = 1 GcV~ 2 , which are smaller, than 
Gcrit- Now, in order to compare our results with those of Wang, it is convenient to divide, as in [l[, the numerical 
calculations into several stages. 



Table I: The case G = 0.5 GeV" 2 . Table H: The case G = 1 GeV~ 2 . 



too [GeV] 





0.0003 


0.003 


0.03 


0.3 




too [GeV] 





0.0002 


0.002 


0.02 


0.2 


to [GeV] 





0.00043 


0.0043 


0.043 


0.4 




to [GeV] 





0.0005 


0.0051 


0.051 


0.4 


2\e\B m [GeV 2 ] 





1.38 


2.01 


3.59 


10.24 




2\e\B m [GeV 2 ] 





0.435 


0.67 


1.36 


4.86 


M{m ,B m ) [GeV] 





0.00194 


0.0154 


0.113 


0.74 




M{m 0l B m ) [GeV] 





0.00212 


0.0165 


0.116 


0.65 


R = M (to , B m )/m 




4.51 


3.58 


2.64 


1.85 




R = M(toq, B rn )/m 




4.16 


3.24 


2.30 


1.64 



(i) First, we put B = and find the dynamical quark mass m, by solving equation (|3|) for different values of the 
bare mass too. For some representative values of too (see the first line in Tables I, II) the corresponding values of m 
are presented in Tables I, II (see the second line there). 

(ii) Next, one should find such a value of the magnetic field B m , for which the solution of equation ([5]) at to = 
coincides with m, i.e. M(m^ — 0, B m ) — to. For the values of to from Tables I, II the corresponding values of B m are 
presented in the third line of Tables I, II. Since A 2 / |e| « 3 GeV 2 « 1.6T0 20 G, we see that in the NJL model the values 



The continuous chiral symmetry remains to be broken due to the presence of the isospin-violating electric charge matrix Q in the 
covariant derivative of the modified Lagrangian (see the remarks above @). 
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of the quantity B m have the order of magnitude of the NJL characteristic magnetic scale B c , B c = A 2 /|e|. In contrast, 
in QED (see [l[) the values of the corresponding quantity B m< , are unrealistically higher than the characteristic QED 
magnetic scale, the Schwinger magnetic field m 2 /|e| « 4.4 • 10 13 G. 

(iii) Finally, for each fixed value of mo and B m we have solved the gap equation ([5]) and found the corresponding 
dynamical quark mass M(mo, B m ) as well as the ratio R = M(rno,B m )/m (these quantities are given in the fourth 
and fifth lines of Tables I, II respectively) which in some sense might serve as a measure of the dynamical quark mass 
enhancement effect 1] . It is seen from these tables that for the considered values of mo we have obtained R < 5 in 
the framework of the NJL model |T]), and even R < 2 for the physically interesting case of a dynamical quark mass 
m = 0.4 GeV (with bare quark mass m = 0.3 GeV for G = 0.5 GeV" 2 or m = 0.2 GeV for G = 1 GeV~ 2 ). For 
comparision, let us quote the value R ss 10, which was obtained in the same manner in QED [1]. 

Judging about the possibility of the enhancement effect in terms of the quantity R only, one might conclude that in 
the framework of the NJL model with a rather weak interaction, G < G cr u , the generation of a dynamical quark mass 
in a strong magnetic field is still enhanced at nonzero bare quark mass (since R ~ 2 in physically reliable cases with 
m w 0.4 GeV). But here this effect is not so pronounced as in QED, where R k, 10. On the other hand, one should 
keep in mind that in the NJL model the enhancement of the dynamical quark mass takes place only at sufficiently high 
magnetic fields B > B m ~ 10 20 G. Indeed, our calculations show that in a more interesting case with realistic values 
of an external magnetic field B < B p h ys = 10 15 G, which are typical values of magnetic fields on the surface of young 
neutron stars, the dynamical quark mass M(mo =/= 0, B) is with great accuracy equal to the dynamical (constituent) 
quark mass m = M(m 7^ 0,0) at B = 0. As a result, we see that for realistic values of B the enhancement effect 
is absent. In contrast, in QED such magnetic fields provide a few percent increase of the dynamical electron mass 
in comparision with m e at B — 0, which is sufficient for observation of this effect in experiments To better 
understand the absence of the enhancement effect in the NJL model at B < B p j lys , one should take into account that 
in the framework of the NJL model the field Bphys is comparatively weak, since Bphys B c , while in QED the field 
Bphys is comparatively strong, since B p h ys S> m 2 /|e|, i.e. it is much greater than the Schwinger field. 

The strong coupling regime (G > G cr n). Since NJL models are considered to be effective theories for low energy 
QCD only at G > G cr u, we have studied the influence of an external magnetic field B on the dynamical quark mass 
also in this case. At G > G cr u and B = 0, the values of the NJL model parameters can be fixed through fitting of 
experimental data, and the typical set of A, mo, G looks like [ll|: A = 0.6 GeV, mo = 0.005 GeV, G — 6.73 GeV -2 , 
which corresponds to G cr u ~ 4.57 GeV -2 and the dynamical quark mass M(mo, B = Q) w 0.4 GeV. 

Now, using the gap equation |5J) , it is possible to conclude that at the value of the external magnetic field B = 
B c = A 2 /|e| w 6.4 • 10 19 G, which is the characteristic magnetic scale of the model for the above chosen value of 
A = 0.6 GeV, the dynamical quark mass M(mo, B c ) exceeds the dynamical quark mass M (mo, B = 0) no more than 
by 20%. At B w 25B C , the corresponding dynamical quark mass is ten times larger than M(mo, B = 0), etc. Hence, 
at G > G cr u, and in a rather strong external magnetic field B > B c the enhancement of the dynamical quark mass 
also takes place. 

However, for values of a magnetic field B smaller than B c , the excess of M(mo,B) over M(nio,B = 0) sharply 
decreases. (Note, that it is just in the region B < B c that the dynamics of QCD is qualitatively similar to that in 
the NJL model [17j].) Indeed, at B = 0.1B C it is equal to 0.3% etc., and for the value B = B p h ys = 10 15 G (the field 
on the surface of young neutron stars) the difference between M(mo, B p h ys ) and M(mo, B = 0) starts from the 10-th 
significant digit, i.e. it is vanishingly small. Therefore, for sufficiently small B < B p h ys <§C B c the enhancement effect 
is absent both in the NJL model and QCD 2 , and hence, in physical applications one might ignore the dependence 
of the dynamical quark mass on an external magnetic field B in this range. In spite of this fact, there are other 
phenomena, which can be observed at B < B p h ys in dense quark matter. Among them are the magnetic oscillation 
effect and other effects [l2|, E3| that are not directly related to the behaviour of the dynamical quark mass vs B. They 
are connected mostly with the thermodynamical properties of the system. 

We remark in conclusion, that in order to answer the question raised at the beginning of the paper, one should 
first establish the ranges for the external magnetic field B under consideration. Then, if B varies in a certain vicinity 
of Bphys = 10 15 G, the considered enhancement effect is intrinsic to QED (since here the magnetic field B p hy S can 
be considered to be strong enough). At the same time, in QCD or the NJL model, the dynamical quark mass is not 
influenced by these realistic values of an external magnetic field. However, if B is rather strong, i.e., B > B c , the 
enhancement of a dynamical fcrmion mass does occur in QED and the NJL-type models, etc. 

Note that in QCD in a strong magnetic field the situation with the enhancement effect might be very involved. 
Indeed, as it was shown in [17] in the chiral limit of QCD, the dynamical quark mass at B > B c behaves quite 
unexpectedly since in a wide range of strong magnetic fields it is suppressed in comparision with the dynamical quark 
mass at B = (one of the reasons is that QCD is an asymptotically free model) . The same behaviour of the dynamical 
quark mass might be inherent to QCD with nonzero bare quark mass. 



2 The behaviour of the quark condensate S(-B) at small values of magnetic fields B <C B c was also considered in the framework of the 
chiral effective theory It is easily seen that in this case at B < B p hys the quark condensate E(-B) exceeds S(0) also in sufficiently 
small fractions of a percent and, hence, the dependence of chiral condensate on B might not be allowed for. 



4 



Acknowledgement 

The authors are grateful to D. Ebert for careful reading of the manuscript and fruitful discussions as well as to V.P. 
Gusynin and S.Y. Wang for useful correspondence and remarks. 



[1] S.Y. Wang, Phys. Rev. D 77, 025031 (2008). 

[2] Yu.M. Loskutov, V.V. Skobelev, Teor. Mat. Fiz. 38, 195 (1979); Vestn. Mosk. Univ. Fiz. Astron. No. 6, 95 (1983); V.P. 
Gusynin and A.V. Smilga, Phys. Lett. B 450, 267 (1999). 

[3] V.A. Miransky, Prog. Theor. Phys. Suppl. 123, 49 ( 1996); A.S. Vshivtsev, B .V. Magnits ky, V.C. Zhukovsky, and K.G Kli- 
menko, Phys. Part. Nucl. 29, 523 (1998); Y.J. Ng, |arXiv:hep-th/9803074l C.N. Leung, |arXiv:hep-th/9806208| K.G. Kli- 
menko, [arXiv:hep-ph/9809218| V.P. Gusynin, Ukr. J. Phys. 45, 603 (2000). 

[4] K.G. Klimenko, Teor. Mat. Fiz. 89, 211 (1991); Teor. Mat. Fiz. 90, 3 (1992); Z. Phys. C 54, 323 (1992); I.V. Krive and 
S.A. Naftulin, Yad. Fiz. 54, 1471 (1991). 

[5] V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, Phys. Rev. Lett. 73, 3499 (1994). 

[6] V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, Phys. Rev. D 52, 4747 (1995); D.K. Hong, Y. Kim, and S.J. Sin, Phys. 
Rev. D 54, 7879 (1996); A.V. Kuznetsov and N.V. Mikheev, Phys. Rev. Lett. 89, 011601 (2002). 

[7] V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, Phys. Rev. Lett. 83, 1291 (1999); Nucl. Phys. B 563, 361 (1999). 

[8] C. N. Leung and S.-Y. Wang, Nucl. Phys. B 747, 266 (2006); Annals Phys. 322, 701 (2007). 

[9] A. Ayala, A. Bashir, A. Raya and E. Rojas, Phys. Rev. D 73, 105009 (2006); Phys. Rev. D 77, 093004 (2008). 
[10] A.S. Vshivtsev, K.G. Klimenko, and B.V. Magnitsky, Nuovo Cim. A 107, 439 (1994); JETP Lett. 62, 283 (1995); Theor. 
Math. Phys. 106, 319 (1996); D.M. Gitman, S.D. Odintsov, and Yu.I. Shil'nov, Phys. Rev. D 54, 2964 (1996); D. Ebert 
and V.C. Zhukovsky, Mod. Phys. Lett. A 12, 2567 (1997); GW. Semenoff, LA. Shovkovy, and L.C.R. Wijewardhana, Mod. 
Phys. Lett. A 13, 1143 (1998); Phys. Rev. D 60, 105024 (2000); E.J. Ferrer and V. de la Incera, Int. J. Mod. Phys. 14, 
3963 (1999); Phys. Lett. B 481, 287 (2000); V.C. Zhukovsky and K.G. Klimenko, Theor. Math. Phys. 134, 254 (2003); 
Theor. Math. Phys. 124, 1132 (2000); R. Gonzalez Felipe, G.M. Marques, and J.E. Ribeiro, Nucl. Phys. A 778, 30 (2006); 
N. Sadooghi and A. Sodeiri Jalili, Phys. Rev. D 76, 065013 (2007); M. d. J. Anguiano-Galicia, A. Bashir, and A. Raya, 
Phys. Rev. D 76 , 127702 (2007); B. Hiller, A. A. Osipov, A.H. Blin, and J. da Providencia, Phys. Lett. B 650, 262 (2007); 
larXiv:0802. 31931 
[11] M. Buballa, Nucl. Phys. A 611, 393 (1996). 

[12] D. Ebert and K.G. Klimenko, Nucl. Phys. A 728, 203 (2003); Phys. Atom. Nucl. 68, 124 (2005). 
[13] V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, Phys. Lett. B 349, 477 (1995). 

[14] D. Ebert, K.G. Klimenko, M.A. Vdovichenko, and A.S. Vshivtsev, Phys. Rev. D 61, 025005 (2000); M.A. Vdovichenko, 

A.S. Vshivtsev, and K.G. Klimenko, Phys. Atom. Nucl. 63, 470 (2000). 
[15] T. Inagaki, D. Kimura, and T. Murata, Prog. Theor. Phys. Suppl. 153, 321 (2004); Prog. Theor. Phys. Ill, 371 (2004). 
[16] H. Bateman and A. Erdeyi. Higher Transcendental Functions. McGrawHill, New York, 1953. 
[17] V.A. Miransky and LA. Shovkovy, Phys. Rev. D 66, 045006 (2002). 

[18] LA. Shushpanov and A.V. Smilga, Phys. Lett. B 402, 351 (1997); N.O. Agasian and LA. Shushpanov, Phys. Lett. B 472, 
143 (200 0); T.D. Cohen, D .A. McGady and E.S. Werbos, Phys. Rev. C 76, 055201 (2007); P. Castelo Ferreira and J. Dias 
de Deus. larXiv:0707.4200l E.S. Werbos, larXiv:0711.2635l 



